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We present a real-time path integral theory for the rate of electron transfer reactions. Using graph theoretic techniques, the 
dynamics is expressed in a formally exact way as a set of integral equations. With a simple approximation for the self-energy, 
the rate can then be computed analytically to all orders in the electronic coupling matrix element. We present results for the 
crossover region between weak (nonadiabatic) and strong (adiabatic) electronic coupling and show that this theory provides a 
rigorous justification for the salient features of the rate expected within conventional electron transfer theory. Nonetheless, we 
find distinct characteristics of quantum behavior even in the strongly adiabatic limit where classical rate theory is conventionally 
thought to be applicable. To our knowledge, this theory is the first systematic dynamical treatment of the full crossover region. 



I. INTRODUCTION 

Since the landmark papers of Marcus QJ^] on the the- 
ory of electron transfer (ET) reactions, the central role of 
the solvent in determining the ET rate has been well rec- 
ognized. The ordinary ET reaction is an activated pro- 
cess , in which there is a free energy barrier separating 
reactants and products. This barrier is the result of the 
often strongly solvated donor and acceptor states, and 
the transfer of an electron then requires a specific large- 
scale reorganization of the environment, usually achieved 
through equilibrium fluctuations. 

ET reactions are commonly thought to be well under- 
stood in two distinct limits (see Fig. |l|). First, for a 
very weak electronic coupling or the so-called nonadia- 
batic limit, the ET can be treated with a perturbation 
theory 0, which to lowest order in the electronic cou- 
pling leads to the Golden Rule expression for the rate. 
For the opposite adiabatic limit in which the electronic 
coupling is large, the rate is thought to be controlled by 
the motion of the solvent on the lower electronic state 
and well described by classical activated rate theory 
In both nonadiabatic and adiabatic theories, there is an 
assumed separation of timescales between the motions of 
the electron and the solvent. In the nonadiabatic limit, 
the solvent fluctuations are assumed fast and the ET is 
controlled by thermal excitations of the solvent to near 
the crossing region followed by a fast transversal of the 
crossing region by the solvent. The ET rate is therefore 
controlled by the probability of an electronic transition 
proportional to the square of the electronic coupling. The 
adiabatic limit on the other hand is characterized by slow 
solvent fluctuations such that the Born-Oppenhcimer ap- 
proximation for the electron can be invoked. Therefore 
according to conventional theories in these two limits, the 
ET rate is expected to depend on the electronic coupling 
7iA/2 like 



, J A 2 cxp(— Eo/kT), nonadiabatic 
ET ~ \ exp(-[£b - hA/2]/kT), adiabatic ' 1 

where T is the temperature and k is the Boltzmann con- 
stant. 

Regarding the assumed separation of timescales, there 
are at least two aspects of ET reactions which re- 
main largely unclear and demand more complete the- 
oretical treatments. First, there is the question of 
the "crossover region" . By this we mean the re- 
gion between the nonadiabatic and the adiabatic lim- 
its in which the electronic matrix element changes from 
the weak coupling to the strong coupling limit. In 
this region, the assumed separation of timescales is 
necessarily violated. In fact, a successful dynamical 
crossover theory that can unify the two limits requires 
inclusion of high order electronic coupling terms |5]-||, 
and such a theory has yet to be found. Second, as 
is true in almost all condensed phase ET reactions, 
the frequency spectrum of solvent motions is charac- 





FIG. 1. Potential surfaces for nonadiabatic (a) and adi- 
abatic (b) electron transfer. In the nonadiabatic case, the 
Landau-Zener mechanism facilitates transfer at the crossing 
point of two diabatic surfaces. The adiabatic transfer is often 
described as a barrier crossing on the lower adiabatic energy 
surface. 
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terized by a broad distribution of timescales. Hence, the 
assumption of timescale separation, while perhaps true 
for parts of the spectrum, may not be so for others. While 
it is commonly believed that fast solvent polarization has 
little effect on the ET rate [||J6), unconventional behav- 
iors have been suggested for solvents with high-frequency 
modes M. Indeed, Stuchebrukhov and Song || have re- 
cently reported that high-frequency components in the 
solvent spectrum can substantially renormalize the in- 
trinsic electronic coupling to slow down the otherwise 
adiabatic rate. The same adiabatic renormalization ef- 
fects have also been recognized by others in a somewhat 
different context ]9|Jl0|. The questions surrounding the 
influence of high-frequency modes in adiabatic reactions 
seem far from settled. 

The search for a unified treatment that can connect 
the conventional nonadiabatic and adiabatic theories 
presents a considerable challenge. Early attempts were 
made by Levich jll] and Zusman to use the Landau- 
Zener formula to incorporate higher-order effects in the 
electronic coupling into the ET rate. Along a similar 
line of thinking, Garg, Onuchic and Ambegaokar used 
an approximate real-time path integral treatment in an 
attempt to account for adiabaticity effects. 

More recent approaches |^,|l4|,|]|ll| have been largely 
based on an imaginary-time formulation of the prob- 
lem. Most noticeable of these is the centroid free energy 
method |u| employed by Gehlen et al. HQ- Exten- 
sions of this to complex-valued centroid coordinates by 
Stuchebrukhov and Song |8| have also been useful for ob- 
taining results for the inverted region where nuclear tun- 
neling is crucial. Another imaginary-time theory which 
does not explicitly employ the centroid formulation has 
been given by Cao et al. |l5| who focus on the compu- 
tational aspects of instanton methods. All these theories 
rely on an assumed relationship between the ET rate and 
an analytic continuation of the partition function for the 
tunneling system |23f| . This relationship, however has not 
been derived from first principles for the ET problem. 

In this paper, we take a rather different approach to the 
problem. With a spin-boson model for the ET reaction, 
our method sets out from an exact path integral formu- 
lation |Io|]l7|| for the real-time dynamics of the electron. 
From the path integrals, it is then easy to see that the ET 
rate is related to clusters of kink/anti-kink pairs (called 
"blips") jl0| which interact with each other through a 
nonlocal influence functional. When the fluctuations of 
the bath polarization are sufficiently fast compared to the 
motions of the electron, these interblip interactions van- 
ish. In this case, the rate reduces to nonadiabatic theory 
and becomes completely equivalent to a low-order per- 
turbation theory in the electronic coupling. On the other 
hand, when the bath motions are sufficiently slow com- 
pared to the electron, interblip interactions force blips 
to bunch up into clusters. In each cluster, multiple elec- 
tronic transitions are possible, and the rate is now given 



by a series containing high-order terms in the electronic 
coupling. Using standard graph theoretic techniques and 
a simple approximation for the self-energy (i.e., the sum 
over all irreducible diagrams), we are able to sum up an 
infinite number of terms in this series and obtain an ex- 
pression for the ET rate not only in the nonadiabatic 
and adiabatic limits, but for the entire crossover region. 
To our knowledge, this is the first successful real-time 
formulation that can span the whole crossover region. 

Section || gives a brief overview of the path-integral 
formulation of the dynamics of the spin-boson system. 
With minor modification, this material follow closely re- 
views by Leggett [|o| and Weiss JTt) . In Section III, 
the formal expression for the time-dependent dynamics 
is transformed diagramma ticall y to yield a rigorous result 
for the ET rate. Sections III A and III B give a pictorial 
interpretation of the results, and Section III C discusses 
general ET behaviors in the two limits i n term s of topo- 
logical features of the graphs. Section [II D| presents a 
simple approximation which allows us to resum the en- 
tire perturbation series for the entire crossover region. 
Some numerical results for the crossover behavior are also 
discussed. Finally, section |y| concludes this paper with 
a brief summary of the most important features of the 
present theory and our results. 



II. DYNAMICS OF ELECTRON TRANSFER IN 
THE SPIN-BOSON MODEL 

For the discussion of ET dynamics in a donor-acceptor 
system in the condensed phase, we take the simple spin- 
boson model which has been used in many previous stud- 
ies. The Hamiltonian is 



where 



H — Hs + Hi + H-g 



H s = --(Aa x +ea z ) 



Hi 



2 ^ 



+ m a uj a x a , 



(2) 

(3) 
(4) 

(5) 



defining two localized electronic states (with a possible 
asymmetry e) and their overlap in the tight-binding ap- 
proximation with a bilinear coupling to a number of sol- 
vent modes, which are assumed harmonic. 

Only two quantities of the oscillator bath are relevant 
to the dynamics of the electron. They are the temper- 
ature associated with the initial thermal density matrix 
of the bath, and its spectral density 



-S(u 



(0) 
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A commonly used form of J(ui), which in the classical 
limit describes a frequency-independent friction, is the 
ohmic spectral density 



— JH = au;e 



(7) 



where a is a dimensionless parameter which character- 
izes the strength of the dissipation. We shall use this 
spectral density throughout this paper. Classically, the 
ohmic bath is characterized by a single parameter - the 
reorganization energy HA — 1cthio c . 

Using standard path-integral techniques [p^|JlC| , [l^ ,|3|| , a 
partial trace can be performed over the oscillator coor- 
dinates to eliminate them from the problem. One is left 
with an effective double path integral expression for the 
dynamics of the reduced density matrix for the electron. 

The forward- and reverse-time paths for the electron 
o~(t) and <r'(t) are piecewise constant functions of time, 
alternating between the two allowed values ±1. This par- 
ticular class of functions can obviously be parametrized 
by the number of transitions between the two allowed val- 
ues, the times at which these transitions occur, and the 
constant values the function assumes in the intervening 
intervals. 

Applying this parametrization, and replacing a(t) and 
o~'(t) with symmetrized and antisymmetrized coordinates 
= {a{t)-a'(t))/2 and n(t) = (<r(f) +er'(f))/2, an ex- 
act path-integral expression for the transition probability 
from position cr, to crj can be written as ||io| , [l7f 

00 ( A 2 

P<Ti<Tf(t) = ti<Ti,<Tf + ^ 

n=l 

* n l-l 

x / v n {tj} n^n^*' ^ 

te,^=±i}/=i fc=i 

where we have used the shorthand notation 

t t tan ta 

J V n {tj} = J dt 2n J dt 2n -i ... J dti (9) 



through 

Bi = exp(-&l-i&lS2i,2l-i + i&i-Mhi - hi-i) 



21,21-1) , 



(11) 



and 



Flk = exp (— ^ 2 l^2k-lS2l,2k-l — ^2l-1^2kS2l-l,2k 
— £,2lS,2kS 2 l,2k —£,2l~l£,2k-lS 2 l-l,2k-l) 
exp (i^2imk-lR2l,2k-l + i&l-lV2k-lR2l-l,2k-l) 
+i&lV2kR2l,2k +i&l-lV2kR2l-l,2k) , (12) 

where S mn = S(t m - t n ) and R mn = R(t m - t n ). Be- 
cause of the two-state nature of the electron coordinate, 
the charges are constrained such that £22-1 = — and 
T]2i = —TI21+1- These constraints imply that there are 
two kinds of time intervals on the double path - those 
intervals with a = a' (dubbed 'sojourns' by Leggett et 
al. 0) and those with a = -a' (called 'blips')- The 
double path alternates between sojourns and blips and 
the expression (||) for P a af {t) involves pair interactions 
between the charges. In (g), we have deliberately broken 
up the charge-charge interactions into "intrablip" terms 
(B) and "interblip" terms (F). 

All of the transition probabilities P CTi(T/ (i) are related 
to each other by the conservation of probability. By sym- 
metry they can be obtained from a single quantity 



P(t) = (a,(t)) = P+,+(t) ~ P+,-(t) 



(13) 



After some (possibly complicated) initial transient, P(t) 
should follow a simple exponential decay with a inverse 
time constant equal to the net reaction rate T. 



III. DYNAMICAL RATE THEORY FOR 
ELECTRON TRANSFER 

In this section, we will discuss the meaning of (|j) us- 
ing diagrams. We begin with the nonadiabatic limit and 
proceed to the more general case at the end of the section. 



for the time-ordered integrals. The {tj} here are the 
times at which transitions are made on the £(i) and r)(t) 
paths and denoting the directions of these transitions by 
£j and rjj (called "charges"), the interaction terms Bi 
and Fik are related to the real and imaginary part of the 
twice-integrated bath correlation function in the follow- 
ing way 



2 

s(t) + iR(t) = !L [du 

nh J 



J(«) 



(1 — cos wt) coth 



■ i sin ut 



Ti(3uo 



(10) 



A. Dynamical rate expression in the nonadiabatic 
limit 



In the limit of nonadiabatic electron transfer, the infi- 
nite series in (||) can readily be resummed using the non- 
interacting blip approximation (NIB A) |l(^,[l7j . Because 
we shall use this later as the basis of an exact expres- 
sion for the reaction rate, we first discuss in detail the 
diagrammatic treatment of the nonadiabatic limit. 

In the NIBA expression P'(t), the factors Fik, which 
denote interactions among charges belonging to different 
blips, are approximated by unity. This turns the multiple 
integrations in (0) into convolution products, which can 
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FIG. 2. Diagrammatic representation of blip-blip interaction terms Gjk for n = 3. Nodes represent blips, straight lines 
sojourns, and arches interaction terms Gjk- The last five diagrams are irreducible, while the first three are not. 

be Laplace transformed. The infinite sum in (||) then = + • 

becomes a geometric series of the Laplace variable A 



1 



P'(A) = T 1 



A-E 



(i) 



(14) 



where 
S«(A) 

4 1} (A) 



' 4 



4 



cir exp [— Ar — S(t) 

+ie(er + 7 ? i?(r))] , 

dr 77 exp [—At — S'(t) 

+^(er + 7?i?(r))] . 



(15) 



(16) 



In the long-time limit, eq. (JIJ) gives an exponential de- 
cay when transformed back. The net reaction rate T is 
obtained as the smallest real solution of the equation 

r + sW(-r) = o. (17) 

This result is closely related to the golden-rule rate 

r F = -EW(A = 0), (18) 

which is just a first approximation to iterative solution 
of @. ' 

Diagrammatically, we can represent this approxima- 
tion by the following series of graphs 



(19) 

The solid line in (|l^) denotes the NIBA transition proba- 
bility P' aa , , which is equal to a sum over diagrams repre- 
senting all possible alternating sequences of sojourns and 
blips. Dashed lines represent sojourns, each of which 
contributes a factor of unity. Dots represent blips, each 
of which contributes a factor of ±(A/2) 2 £?/, correspond- 
ing to two electronic transitions and the interactions of 
the associated charges. For each internal sojourn line 
the summation over its 77 label is implied. Equally, each 
node also represents a summation over the blip's £ label 
and integration over the two transition times. Switching 
to the Laplace transform eliminates one of these integra- 
tions, and changes the other into an integration over the 
blip time fil| 77 = t% — t%-\. These integrals factorize, 
allowing the resummation (|l4|). 



B. Exact dynamical rate expression 

We now turn to an exact diagrammatic representation 
for Pit). Using NIBA as the zero-order solution, we de- 
velop a perturbation series for P(t). In place of the NIBA 
Fik = 1, employing the ansatz 



4 



1 



Ik 



1 + G 



Ik 



(20) 



expands the product Y[ Flk m (@) into a sum, which is 
conveniently represented by another series of diagrams. 
To generate all terms in this series, we take (|l^) and add 
to it all diagrams in which each pair [l, k) of nodes I ^ k 
may now be connected by (at most) one curved line that 
represents the interaction Gjk- As an example, all third 
order diagrams generated this way are depicted in Fig. [|. 

Some of these diagrams are reducible in the sense that 
they can be separated into disjoint parts by removing a 
sojourn line. (An isolated single node is taken as the 
smallest connected diagram). Every reducible diagram 
can thus be uniquely decomposed into a sequence of ir- 
reducible diagrams and sojourn lines. The linked- cluster 
theorem similarly states that the sum of all reducible di- 
agrams is graphically represented by 



-O 



-0---0- 



(21) 

where the linked-cluster sum (or self-energy, in the lan- 
guage of quantum field theory) represented by the open 
circle is 



o 






(22) 



In (pi]), the double line represents the exact P(t) and 
dashed lines are again sojourns. The smallest irreducible 



diagram in fl22j ) is just a single node (i.e. the right-hand 
side of the first line of (|2|) , representing a blip that does 
not interact with anything else. 

In this decomposition each of the multiple integral in 
(ph in the Laplace transform p(X) becomes a convolution 
product of several terms, each corresponding to an ir- 
reducible self-energy diagram. Formally, applying the 
linked-cluster theorem means rearranging the summa- 



tions in (g) and (20) to make full use of this property 
of the integrand. The new outermost summation is then 
taken over the number of connected subdiagrams and 
leads again to a geometric series for p(X), 



i 



(23) 



The self-energies 

E,(A) = \ J2 ^no , £„(A) = \ E »» E * ( 24 ) 

i7o = ±l »jo = ±l 

are the symmetric and antisymmetric part of the sum 
over all irreducible diagrams, and 



00 / A 2 \™ r 

n-1 



1 E - 



Bie~ XT1 



(25) 



\[< v - E II 

j=i cec„ (j.k)ec 

Here C n denotes the set of all irreducible diagrams with n 
nodes, and the notation (j, k) 6 C means that the graph C 
contains a line connecting nodes j and k. The integration 
variables have been changed to the blip times Tj = t%j — 
t-zj-i and sojourn times Sj = t2j+i — t%j, where t\ now 
marks the beginning of an irreducible cluster. 

In this expression, an irreducible diagram may contain 
any number of internal nodes that are not part of any 
interaction line. Such subdiagrams enclosed by two 'in- 
teracting' nodes may be summed over in a fashion similar 
to NIB A and treated as an 'extended sojourn'. Equation 
(S) then assumes the equivalent form [M 



o 



+ 






(26) 



5 



To obtain this series, we eliminate from ( j22] ) all contigu- 
ous sequences of unconnected internal blips and replace 
each sequence by a single NIBA (solid) line. For exam- 
ple, the first line in (p6| ) in the result of resumming the 
first three lines in (|22|). Mathematically, this translates 
to 



E 

CGG 



A 2 

' 4 



n(C) 



dT n l ds n - 1 




n(C)-l 

e ^- ati n e 

{Sj, % =±i} i=i ^-=±1 

xe -A(, 3 +r 3+1 ) jB / +i -Q G . fc) 

y,fc)ec 



(27) 



where n(C) is the number of nodes in a graph C, and 
P^. cr/ (t) is the NIBA transition probability. The terms 
B'j differ from Bj by the substitution rjj — > rjj. It is im- 
portant to note that the set G contains only irreducible 
diagrams with a further restriction: Except in the trivial 
single-blip diagram, each node must be part of an inter- 
action line |2~i| . 

The reaction rate T can now be obtained from either 
( p5[ ) or (|27j ) as the solution of 



r + s s (-r) = o 



(28) 



For relaxation rates slow compared to the dynamical 
timescales of the system, (|2^) reduces to 

(29) 



S S (A = Q). 



With (^5|) and ( p7j ) we have thus found a rigorous re- 
sult for the dynamics of the spin-boson system. Together 
with reasonable approximations for the self-energy, we 
can produce systematic improvements on the NIBA. 



C. Leading adiabatic corrections: breakdown of 
steepest descent 

Armed with these formulae, we can study the adia- 
batic corrections, i.e. deviations from the 0(A 2 ) NIBA. 
A startling feature of the adiabatic electron transfer prob- 
lem is the fact that the Gaussian approximation for the 
integrand in (|25|) obtained by a high-temperature expan- 
sion of the correlation functions S and R breaks down 
for terms of any order higher than A 2 . A simple exami- 
nation of the fourth-order term will reveal the problem. 
Assuming for the sake of argument that temperatures 
in the classical regime kT 3> Tiuj c justify an expansion of 
the functions S(t) and R(t), the fourth-order term of 
takes the form 



,(2) _ 
J V0 - 



A 4 



E 



drx I ds\ I dr 2 e 



-Ati — Asi — At2 



x exp(-AfcTT 1 2 /?i + i(e + no^m + Hi) 
x exp(-AfcTT 2 2 /?i + i(e + r7iA)£ 2 T 2 + 1Z 2 ) 
x [exp{-2^ 2 AkTnT 2 /h + i(m - t?i)6At 2 + ft z ) 
-1] . (30) 

The terms 1Z\ , 1Z 2 , and IZi denote contributions of higher 
than quadratic order in the exponent, and h~A = 2ah~uj c 
is the reorganization energy of the environment. With 
Tj = £jTj and Sj — rjjSj, the integrations in 

oo , oo oo 

Z$=A 4 J d~ Sl e-^l J dhe^ J df 2 e-^ 

— -oo — oo — oo 

x exp [AAT(n + f 2 ) 2 /h + % A + e)(fi + f 2 ) + V] 

-EgM 1 )} (31) 



become formally Gaussian for A, TZ — > 0. Due to the de- 
generacy of the quadratic form in the exponent, however, 
one is forced to retain the higher-order corrections TZ in 
order to prevent the integral from diverging. The same 
picture, with the quadratic part of the exponent depend- 
ing only on one single variable r* = £jtj and 2n — 2 
'false' zero modes, emerges for terms of arbitrary order 
A 2 ™. Thus, it is not surprising that a numerical eval- 
uation of the complete expression (|3l| ) leads to results 
quite different from those obtained when approximating 
the bath correlation function by its short-time behavior. 

Following Garg et al. || , a convenient quantitative de- 
scription of the transition from nonadiabatic to adiabatic 
behavior can be given in terms of an adiabaticity factor 
g(A, e, u) c , a, T) defined by 



r-r B 



9 = 



(32) 



A=0 



which specifies the relative weight of the non-trivial di- 
agrams in (p2|), i.e., the multiblip contributions to J25|) . 
When adiabaticity effects reduce the rate, the sum asso- 
ciated with these higher-order diagrams can exceed the 
rate itself in absolute value, making g a positive number. 
A negative g, on the other hand, would indicate higher 
order terms promoting electron transfer. 

Restricting (^5|) to terms of order A 4 for small adia- 
batic corrections g <C 1, (p3) turns into 



.9=1 



r 

1 or 



and g scales as 



A 2 

g = -^g{a, kT/huj c , e/w c ) ■ 



(33) 



(34) 



We have evaluated g(a, kT/?uv c ) numerically for a sym- 
metric system, with results shown in Fig. [| For high 



G 



temperatures kT ^> fiA the adiabaticity factor is only 
weakly dependent on a, but still varies significantly with 
temperature. Here we are in disagreement with Ref. 
Q , where a temperature-independent adiabaticity factor 
g cx a^ 1 was found 

At lower temperatures, but still with kT > huj c we find 
that the function g(a, kT/huj c ) obeys a further approxi- 
mate scaling relation 



.</(". kT/Ji^,.) \ ( -j^ 



hA 

kT 



(35) 



where k ks 0.8. The function w(x) approximately follows 
an Arrhcnius law with a very small activation energy 
w 0.065 fiA. Fig. ^ shows numerical data for the universal 
function w(x) with a varying from 4 to 100, and with 
temperatures from kT /Tuv c — 4 to 80. 

Our result ( |35| ) also shows a significant discrepancy 
compared to thermodynamic rate expressions as, e.g., 
in the approach by Song and Stuchebrukhov || using 
Langer's method p3[] , where the only relevant parameter 
for an environment of 'classical' modes with Tiio c -C kT 
are the reorganization energy TiA and temperature T. In 
their approach, the characteristic frequency lu c ceases to 
be a relevant parameter. The scenario most frequently 
invoked to justify Langer's method, however, supposes 
that the action can be factorized near a saddlepoint, with 
a unique unstable coordinate being the reaction coordi- 
nate, and other coordinates undergoing only fluctuations. 
To our knowledge, this property of the action has not 
been demonstrated for adiabatic ET or for the crossover 
region, nor has the unstable coordinate been identified. 



D. Theory for the crossover from nonadiabatic to 
adiabatic electron transfer 



Although approximating (25) by a finite number of di- 
agrams has provided us with valuable information about 
the leading corrections to nonadiabatic theory, a better 
approximation is needed for the crossover to adiabatic 
electron transfer. In the delocalized transition state, the 
electron may oscillate back and forth many times before 
settling on one of the two sites. Moreover, in the case of 
strong dissipation correlated recrossings of the reaction 
barrier must be considered for an accurate rate expres- 
sion. 



An exact resummation of (25) or (27) seems out of 
reach without further insight into the formal structure of 
the problem. We shall present an approximate resumma- 
tion of ( p7j ) instead, which captures the most important 
contributions in several parameter regions. Our approx- 
imation is guided by the fact that the interaction terms 
Gjk fall off rapidly when the separation of blips is larger 
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FIG. 3. Adiabaticity factor g obtained from the lead- 
ing corrections to the Golden Rule. Temperature values are 
kT/hw c — 4 (solid line), 6 (long dash) , 8 (short dash), and 
10 (dotted line). The adiabaticity factor of Garg et al. 
which is temperature-independent, is given for comparison 
(dash-dotted line) . 
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which is typically a much shorter timcscalc 



than the inverse of the reaction rate. For any diagram 



10 15 20 

A/kT 



FIG. 4. Numerical data for the universal function 
w(hA/kT). Temperature varies from kT/hui c = 4 to 80, with 
a ranging from 4 to 100. Dotted line: fit by a exp(— b HA/kT) 
with a = 0.63 and b — 0.065. The dashed line interpolates 
between data points. 

with nested interaction lines, this cutoff applies to the 
sum of all blip and sojourn times spanned by the outer- 
most interaction line. This constraint decreases the vol- 
ume of the integration domain by a factor n\ compared to 
case of independent cutoffs for each of n integration vari- 
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ables. Diagrams with nested lines will therefore not be 
taken into account, and diagrams with crossed lines will 
similarly be neglected. Therefore, in the self-energy we 
only retain a set of 'bridge diagrams', which correspond 
to only the first two lines of (p6|). 

It is now useful to keep the integration over the last 
blip time r and the summation over its indices £ and rj 
separate from the preceding blips. The self-energy terms 
can then be written as a power series of the form 



oo 

S„ =2£ dr A 2 "B n (r,£,%) 

*, — 1 J f l_1 



(36) 



n=ln f=±l 



For the bridge diagrams, the terms B n obey a linear re- 
cursion relation 

oo 



xB„(r',e',7/) , (37) 



where 



Bi(r,tv) = - i e- AT exph5(r) + ze(er+ryi?(T))] . (38) 
The kernel of the integral transform in ( |37| ) is given by 

oo 

K(T,£,ri;T',e,r l ') = B 1 (T,Z,ri) J dse~ Xs 

o 

x E W'^W^r,^',^'^";*) (39) 
n"=±i 

where 

G(t, t', £ rf, rf\ s) = exp[^^(i?(r'+s +r) - R(t'+ s))} 
x exp[-ir ? "^(i?(r' + s) - i?(s))] 
x exp[-^'(^(r' + s + r) + 5(s) 

-S(s + t)-S(t' + s))] - 1 . (40) 

^From this recursion relation an integral equation for the 
entire series 



(41) 



n=l 



is easily derived, 

r(T,€ >)7 )-A 2 fli(T,€,»j) 



= A 2 f dr' £ *"(t, £, »?; t', rf)Y(/, g, rf) . (42) 

Finally, the linkcd-cluster sum T, Vo is obtained by apply- 
ing the final integration and summation 



oo 

2 J dr Yl Y ( T ^^0) ■ 



(43) 



With (||), @, and @ we are now left with a finite 
number of quadratures and an integral equation, all of 
which can be solved with very modest numerical cost. 

A certain degree of control over our approximation of 
the linked-cluster sum can also be exercised by the inclu- 
sion of a large class of nested diagrams. This can be done 
by replacing the NIBA propagating function i^y(sj-) 

in ( p6| ) by the full expression P^^'^Sj) an d solving the 
equation system given by and (||^) self-consistently 
by iteration. The results we present in this paper do 
not change when these additional diagrams are included. 
An appealing aspect of this approach compared to other 
dynamical methods is the fact that the dynamics at arbi- 
trarily long times can be studied by evaluating the linked- 
cluster sum at small A. 

First, in the high-temperature (almost barrierless) 
case, the ET rate should depend on the electron cou- 
pling A/2 only through the adiabaticity prefactor. In 
this case, electron transfer rates obtained by the dia- 
grammatic method compare very favorably with recently 
published quantum Monte Carlo data [^4|, as shown in 
Fig. ||, which shows the dependence of the rate on the 
electronic coupling. For the parameters we have chosen 
here, the activation factor is close to unity, i.e., the ob- 
served variation of the rate is mostly due to dynamical 
effects. Our results are in agreement with the traditional 




2.0 3.0 

A/co c 

FIG. 5. Comparison of the present theory and QMC simu- 
lations |24| in the crossover region between nonadiabatic and 
adiabatic electron transfer in the high-temperature (almost 
barrierless) case for A./uj c = 4 and kT /Tiuj c = 4. Dashed line: 
Golden Rule expectation. Solid line: Summation over bridge 
diagrams. Symbols: Quantum Monte Carlo data. 
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15.0 




#co c I kT 



FIG. 6. Reactant decay for a symmetric self-exchange 
ET reaction (zero long-time limit) and for an asymmetric 
donor- accept or pair (e = 4Wc, negative long-time limit) for 
the same parameters in Fig. fel Solid lines represent the re- 
summation of the bridge diagrams in ^); dashed lines give 
the noninteracting-blip approximation for comparison. 

expectation of an initially quadratic dependence on the 
electronic coupling that rises and saturates to a plateau 
value for very large electronic couplings. 

For these parameters, as in MC simulations, the time- 
dependent decay of Pit) shows initial oscillations before 
approaching a limit of exponential decay, Fig. [| Note 
that our method correctly reproduces the equilibrium 
value 



Foe. = lim Pit) = tanh 



he 

2kf ' 



(44) 



i.e., the principle of detailed balance is preserved. 

Next, we examine the diagrammatic results in the low- 
temperature (activated) region. In Fig. [?] we present a set 
of Arrhenius plots for increasing values of the electronic 
coupling, but keeping the classical reorganization energy 
constant. At the smallest value A/uj c = 1, deviations 
from the NIBA are minute. For somewhat larger values 
of A, however, one finds significant discrepancies due to 
adiabaticity effects. 

The detailed features of these deviations conform to a 
number of known or intuitively expected characteristics 
of adiabatic electron transfer, (a) The slope of each Ar- 
rhenius plot decreases with increasing A. This of course 
reflects the lowering of the activation free energy through 
the electronic coupling p3|^4T Jq]. (b) For high enough 
temperatures, the absolute value of the rate is signifi- 
cantly reduced from the Golden Rule rate, as predicted 
by dynamical theories 12 l|. (c) All results show some 
degree of nuclear tunneling, visible in a deviation from 



FIG. 7. Rates obtained in the low-temperature (activated) 
case for A/uj c = 30. Results from summation over bridge 
diagrams for A/cu c — 1 (squares), A/ui c — 2 (diamonds), and 
A/u) c — 3 (triangles). The dashed lines indicate the Golden 
Rule rate for comparison (A/o; c = 1, 2, and 3, bottom to 
top). The classical nonadiabatic rate for A/ui c = 1 is given 
for comparison (dotted line). At high enough temperatures, 
Arrhenius-like behavior is observed, as indicated by the match 
between data and straight solid lines. 

Arrhenius behavior at low temperatures, with a rate that 
becomes much higher than its expected classical value. 

We can put observations (a) and (b) on a more quanti- 
tative basis by separating the rate into an Arrhenius-like 
part and a temperature independent prefactor as 



T = vexpi-F*/kT) 



(45) 



Using (45) as an ansatz to fit our data between kT j — 
0.6 hio c and kT I — huj c , we can extract the A-dependence 
of v and F* separately, as shown in Fig. ||. The adia- 
baticity prefactor v shows the expected transition from 
quadratic to saturation behavior, and the activation free 
energy F* shows approximately the expected linear de- 
cline. 

Our diagrammatic theory reproduces the essential 
tenets of both the thermodynamic approach — the mod- 
ification the reaction barrier through the electronic 
coupling — and of dynamical theories — the slowing of the 
reaction due to recrossing effects — in a unified descrip- 
tion. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have presented a novel dynamical ap- 
proach to the theory of electron transfer reactions. The 
spin-boson model is used to represent electron transfer in 
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0.0 1.0 2.0 3.0 0.0 1.0 2.0 3.0 

A / co c A / co c 

FIG. 8. (a) Activation free energy F* as a function of electronic overlap A, using data from Fig. |?] between kT = 0.6 hu> c 
and kT = hcj c . The open square at A = denotes the Golden Rule value, (b) Adiabaticity prefactor v from the same data. 



condensed-phase from a donor to an acceptor site. Re- 
lying on the established real-time path-integral represen- 
tation of the spin-boson dynamics, we have for the first 
time rigorously deduced an exact rate expression in terms 
of the spin-boson parameters. In the nonadiabatic limit, 
our result reduces to the simpler NIBA description of 
spin-boson dynamics. Using the linked-cluster theorem, 
we have found a systematic way of summing up all cor- 
rections to the NIBA rate expression. Since we are using 
a dynamical treatment directly, our method is valid un- 
der rather weak assumptions and it is applicable in any 
parameter region that would yield a long-time limit of 
exponential decay of the reactant population. 

We have demonstrated the relevance of our central re- 
sults and the continued need for systematic groundwork 
for the electron transfer problem by pointing out signif- 
icant discrepancies between results based on the linked- 
cluster sum ( p5| ) and previous approximate real-time or 
imaginary-time approaches. We have also shown the use- 
fulness of our approach by presenting a resummation of 
the most important terms in the linked-cluster sum (27). 
This resummation imposes no approximation on terms 
up to order A 4 , and it retains a large, dominant class of 
diagrams for any order of the electronic or environmen- 
tal coupling. Our results compare very favorably with 
published quantum Monte Carlo data and the results of 
a refined version of the resummation method. 

The present approach provides a unified treatment of 
several characteristic features of adiabatic electron trans- 
fer which were reported earlier, but separately, using in- 
compatible methods. Our results show two competing 
adiabatic effects. First, there is a rate reduction (com- 
pared to the nonadiabatic case), which has previously 



been described in terms of correlated recrossings of the 
Landau-Zener region. This effect is dominant at high 
temperatures. For low enough temperatures, one sees 
a transfer rate that is enhanced over the nonadiabatic 
result. This reflects the lowering of the reaction bar- 
rier by the electronic coupling and by nuclear tunneling. 
When higher demands on computational resources can 
be tolerated, our recursive resummation method may be 
extended to include larger classes of diagrams if higher 
accuracy is desired. 
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